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^ ■ ABSTRACT 

Qh, The recent discovery of a 2Mq binary millisecond pulsar (jPemorest et al.ll2010l ) has 

not only important consequences for the equation-of-state of nuclear matter at high 
Jl^ ' densities but also raises the interesting question if the neutron star PSR J1614— 2230 

^ , was born massive. The answer is vital for understanding neutron star formation in 

^ ' core collapse supernovae. Furthermore, this system raises interesting issues about the 

nature of the progenitor binary and how it evolved during its mass exchanging X-ray 
phase. In this paper we discuss the progenitor evolution of PSR J1614— 2230. We have 
performed detailed stellar evolution modelling of intermediate-mass X-ray binaries 
\^ \ undergoing Case A Roche-lobe overflow (RLO) and applied an analytic parameteriza- 

■ tion for calculating the outcome of either a common envelope evolution or the highly 
I super-Eddington isotropic re-emission mode. We find two viable possibilities for the 
, formation of the PSR J1614— 2230 system: either it contained a 2.2 — 2.6 Mq giant 
• donor star and evolved through a common envelope and spiral-in phase or, more likely, 
I it descended from a close binary system with a 4.0 — 5.0 Mq main sequence donor star 

via Case A RLO. We conclude that the neutron star must have been born with a mass 
of ^1.95 Mq or 1.7 ± 0.15 Mq, respectively - which significantly exceeds neutron star 
birth masses in previously discovered radio pulsar systems. Based on the expected 
neutron star birth masses from considerations of stellar evolution and explosion mod- 
' els, we find that the progenitor star of PSR J1614— 2230 is likely to have been more 

^ , massive than 20 Mq. 
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1 INTRODUCTION 

Neutron stars are formed as compact remnants of massive 
stars (10 — 30 Mq ) which explode in supernovae at the en d of 
their stellar life l|Wooslev et al ] |2002l : lHeger et al.ll2003l '). In 
order to better understand the mechanisms of the electron 
capture and core collapse supernovae knowledge of the dis- 
tribution of birth masses of neutron stars is vital. However, 
in order to weigh a neutron star it must be a member of a bi- 
nary system. This introduces an uncertainty in determining 
the original birth mass of the neutron star since these neu- 
tron stars are often observed in X-ray binaries or, at a later 
stage, as recycled pulsars and hence after they have under- 
gone a phase of mass accretion from their companion star. 
The most precisely measured masses of neutron stars are ob- 



tained in double neutron star systems via general relativistic 
effects. The related post-Keplerian parameters include peri- 
astron advance, redshift/time dilation, o rbital period deriva- 
tive and Shapiro delay (e. g. IWilll |200£ ) . Sha piro delays of 
radio signals from pulsars ( Stairs et al. 11998 ?) have the ad- 
vantage of being measurable also in low eccentricity systems 
if the orbital inclination is such that the pulses passes in the 
vicinity of its companion. This method yields the opportu- 
nity to weigh both neutron stars accurately - and hence also 
determine the mass of the last formed neutron star which has 
not accreted any material. So far, such measurements have 
revealed that even the most massive of these neutron stars 
the non-re cycled pulsars) do not exceed a mass of 1.39 
Thorsett fc Chakrabartv.. 1999. : .Schwab et al.,.20ia ). There 
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is, however, some evidence from neutron stars in X-ray bina- 
ries, e.g. Vela X-1, that suggests neutron stars can be born 
more massive than this value. 

Binary millisecond pulsars are known to be key sources 
of research in fundamental physics. They host the densest 
matter in the observable Universe and possess very rapid 
spins as well as relativistic magnetospheres with outflowing 
plasma winds. Being ultra stable clocks they also allow for 
unprecedent ed tests of gravitation al theories in the strong- 
field regime ijKramer fc Wexl2009h . Equally important, how- 
ever, binary millisecond pulsars represent the end point of 
stellar evolution, and their observed orbital and stellar prop- 
erties are fossil records of their evolutionary history. Thus 
one can use binary pulsar systems as key probes of stellar 
astrophysics. 

Recent Shapiro delay measurements of PSR J1614— 2230 
l|Demorest et al.l [iolO) allowed a precise mass determina- 
tion of this record high-mass pulsar (neutron star) and 
its white dwarf companion. Characteristic parameters of 
the system are shown in Table [T] It is well established 
that the neutron star in binary millisecond pulsar systems 
forms first, descending from the initially more massive of 
the two binary stellar components. The neutron star is 
subsequently spun-up to a high spin frequency via accre- 
tion of mas s and angular momentum once the secondary 
star e volves (Alpar ct al.'l982':'Radhakris hnan &: SrinivasanI 
1 19821 : iBhattacharva fc van den Hcuvcl J^9l|). In this recy- 
cling pha se the system is observable as a low-mass X-ray bi- 
nary (e.g. lNagase|[l989l ) and towards the end of this phase as 
an X-ray millisecond p ulsar (jWiinands &: van der KMI 19981 : 
[Archibald et al.ll2009l ). Although this formation scenario is 
now commonly accepted many aspects of the mass-transfer 
process and the accretion physics (e.g. the accretion effi- 
ciency and the details of non-conservative evolutio n) are still 
not well understood (|Lewin fc van der KlidlioO^ ). 

In this paper we investigate the progenitor evolution of 
PSR J1614— 2230. We are mainly focusing on the important 
X-ray binary phase starting from the point where the neu- 
tron star has already formed. However, we shall also briefly 
outline the previous evolution from the zero-age main se- 
quence (ZAMS) binary until this stage since this evolution 
is important for the birth mass of the neutron star. In Sec- 
tion [2] we discuss the three different possibilities for mass 
transfer toward a neutron star from an intermediate-mass 
star of 2.2 - 5.0 M©, for the Roche-lobe overflow (RLO) 
Cases A, B and C. The evolution of the original ZAMS bi- 
nary until the X-ray phase is briefly discussed in Section [S] 
In Section 0] we compare our resul t s wit h the outcome of 
the independent work bv lLin et al.l (|201lf ) and also discuss 
our results in a broader context in relation to neutron star 
birth masses predicted by stellar evolution and supernova 
explosions. Our conclusions are given in Section [S] 
In Paper II l| Tauris et al.|[201ll ) we continue the discussion 
of PSR J1614— 2230 in view of general aspects of accretion 
onto neutron stars during the recycling process of millisec- 
ond pulsars. 



2 MASS TRANSFER IN X-RAY BINARIES 

Consider a close interacting binary system which consists of 
a non-degenerate donor star and a compact object, in our 



Table 1. Physical parameters of t he binary m illisecond pulsar 
PSR J1614-2230 (data taken from (Demorest ct al. 2010). 



Parameter 



value 



Pulsar mass 
White dwarf mass 
Orbital period 

Projected pulsar semimajor axis 
Orbital eccentricity 
Inclination angle 
Dispersion-derived distance 
Pulsar spin period 
Period derivative 



1.97 ± 0.04 M0 
0.500 ±0.006 M0 
8.6866194196(2) days 
11.2911975 light sec 
1.30 ±0.04 X 10-** 
89.17 ±0.02 deg. 

1.2 kpc 
3.1508076534271 ms 
9.6216 X 10^21 



case a neutron star. If the orbital separation is small enough 
the (evolved) non-degenerate star flUs its inner common 
equipotential surface (Roche-lobe) and becomes a donor star 
for a subsequent epoch of mass transfer toward the, now, 
accreting neutron star. In this phase the system is observed 
as an X-ray binary. When the donor star flUs its Roche- 
lobe it is perturbed by removal of mass and it falls out of 
hydrostatic and thermal equilibrium. In the process of re- 
establishing equilibrium the star will either grow or shrink 
- depending on the properties of its envelope layers as dis- 
cussed below - first on a dynamical (adiabatic) timescale 
and subsequently on a slower thermal timescale. However, 
any exchange and loss of mass in such an X-ray binary sys- 
tem will also lead to alterations of the orbital dynamics, via 
modifications in the orbital angular momentum, and hence 
changes in the size of the critical Roche-lobe radius of the 
donor star. The stability of the mass-transfer process there- 
fore depends on how these two radii evolve (i.e. the radius 
of the star and the Roche- lobe radius). The various possible 
modes of mass exchange and loss include, for example, direct 
fast wind mass loss, Roche-lobe overflow, with or without 
isotropic re-emission, and common envelope evolution (e.g. 
Ivan den Heuvellll994l ; ISoberman et al. I ll997l . and references 
therein). The RLO mass transfer can be initiated while the 
donor star is still on the main sequence (Case A RLO), dur- 
ing hydrogen shell burning (Case B RLO) or during helium 
shell burning (Case C RLO). The corresponding evolution- 
ary timescales for these different cases will in general proceed 
on a nuclear, thermal or dynamical timescale, respectively, 
or a combination thereof. We now investigate each of these 
three cases with the aim of reproducing the parameters of 
PSR J1614-2230. 



2.1 



Case C RLO 
transfer 



dynamical unstable mass 



10^ 



Donor stars in systems with wide orbits (Porb 
10'^ days) prior to the mass-transfer phase develop a deep 
convective envelope as they become giant stars before fill- 
ing their Roche-lobe. The response to mass loss for these 
stars with outer layers of constant low entropy and negative 
adiabatic mass-radius exponents (^ = d\n R/ d\n M < 0) is 
therefore expansion which causes the stars to overfill their 
Roche-lobes even more. To exacerbate this problem, binaries 
also shrink in size if mass transfer occurs from a donor star 
somewhat more massive than the accreting neutron star. 
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Figure 1. Stellar core mass as a function of the initial ZAMS mass, calculated at different evolutionary epochs. Mass transfer initiated 
either before the TAMS, the tip of the RGB or the AGB, corresponds to RLO Case A, Case B and Case C, respectively. Filled circles 
represent models without convective core-overshooting. The open circles show a few examples of core masses at the tip of the AGB 
assuming a convective core-overshooting parameter of <5ov = 0.20. The horizontal red lines indicate the measured mass interval (within 
3-(T error bars) of the white dwarf in PSR J1614— 2230. This white dwarf descends from the core of the donor star in the X-ray binary. 



This causes further overfilling of the donor star Roche-lobe 
resulting in enhanced mass loss etc. This situation is clearly 
a vicious circle that leads to a runaway mass transfer and 
the formation of a contact binary with a common envelop e 
(CE) followed by a spiral- in phase, e.g. iPaczvnskil l|l976f ). 
Ilben fc Liviol (|l993l '). 

A simple estimate of the reduction of the orbit can 
be found by equating the binding energy of the envelope 
of the AGB giant donor to the difference in orbital energy 
before and after the CE-phase. The idea is that the motion 
of the neutron star, once captured in the CE, results 
in friction and thus dissipation of orbital energy which 
can be u s ed to expel the CE. Follow ing the formalism of 
IWebbinkI l|l984l ') and Ide Kooll (|l99(f), the binding energy 
of the envelope at the onset of RLO mass transfer can be 
written as: — G'M2Monv/(A -R2), where M2 is the mass of 
the donor star, with envelope mass Monv, and R2 ~ Rl 
is the Roche-lobe radius of the donor star at the onset of 
the CE-phase. This radius is often calcul ated in terms o f 
its dimensionless Roche-lobe radius, l|Eggletonl Il983l ) 
such that R2 — Rh = ao ■ tl, where ao is the initial orbital 
separation. 

The total binding energy of the envelope includes both 
the negative gravitational binding energy and the positive 
thermal energy. Besides from the thermal energy of a sim- 
ple perfect gas, the latter term also includes the energy of 



radiation, terms due to ionization of H and He and dissoci- 
ation of H2, as well as the cont ribution from the Ferm i en- 
ergy of the degenerate electrons (|Han et al.lll994l . [l995l ). The 
value of the A-par ameter can t hus be calculated from stellar 
structure models (jPewi fc Tau ris 2000, 2001; Tauris & DcwJ 
I2OOII : IXu fc Ll2010bl lal: lLoveridge et al.ll20ia : IIvanoval201ll ') . 
Given the radius of the donor star and the A-parameter 
enables one to estimate the change in orbital separation 
as a result of the neutron star spiral-in and ejection of 
the envelope. Let rice describe the efficiency of ejecting the 
envelope via drag forces, i.e. of converting orbital energy 
{Eorb = — GM2MNs/2a) into the kinetic energy that pro- 
vides the outward motion of the envelope: -Ecnv = ?7cg Ai?orb 
and one finds the well-known expression for the ratio of the 
change in orbital separation: 

a_ _ McorcMNS 1 QN 

ao ^ M2 MnS + 2Menv/(7?ceArL) 

where Mcorc = M2 — Menv is the core mass of the evolved 
donor star (essentially the mass of the white dwarf to be 
formed); Mns is the mass of the neutron star and a is the fi- 
nal orbital separation after the CE-phase. Strictly speaking, 
when considering the energy budget the "effective efficiency 
parameter" should also include the excess energy of the 
ejected matter at in finity - although this effect is probably 
small. Recent work l|Zorotovic et al.l I2OI0I : Ide Marco et al.l 
l201lh suggests that the efficiency parameter is of the order 
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30%, i.e. TTjcc — 0.3, although its uncertainty is large. The 
value may not be universal and could, for example, depend 
on the stellar mass ratio in a given binary. 
During the very short spiral-in and ejection phase of a com- 
mon envelope evolution (~ 10^ yr) it is a good approximation 
to assume that the neutron star does not accrete any sig- 
nificant amount of matter given that its accretion is limited 
by the Eddington luminosity corresponding to a maximum 
accretion rate of lO~*M0yr~^, depending on the exact 
chemical composition of the accreted material and the ge- 
ometry of its flow. 

The possibility of hypercritical accretion onto the neu- 
tron star during the spiral-in phase has been suggested to 
lead to significant mass increas e and possible collapse of the 
neutron star into a black hole (jChevaheijllQoi : iBrownll 19951 : 
iFrver et al.lll996l ). However, there is observational evidence 
that, at least in some cases, this is not the case. The recent 
determination of the low pulsar mass in PSR J1802— 2124 of 
1.24 ± 0.11 Mq (|Ferdman et al.ll2010l ) clearly demonstrates 
that this 12.6 millisecond recycled pulsar did not accrete 
any significant amount of matter. This pulsar is in a tight 
binary with an orbital period of only 16.8 hours and it has 
a carbon-oxygen white dwarf (CO WD) companion of mass 
0.78 ± 0.04 Mq. With such a small orbital period combined 
with a massive white dwarf there seems to be no doubt that 
this system evolved through a common envelope and spiral- 
in phase and the low pulsar mass reveals that very little 
mass has been accumulated by the neutron star during this 
phase. However, the recycling of this pulsar does require 
some 10~^ Mq of accreted material (see further discussion 
in Paper II). While in principle, PSR J1802-2124 could have 
formed via accretion-induced collapse of an O-Ne-Mg WD, 
so far there are no models which suggest this. 

Before we proceed to discuss the case of 
PSR J1614— 2230 let us introduce a new parameteri- 
zation for calculating the outcome of a common envelope 
evolution. We define a mass ratio parameter k = qo/q where 
go and q represent the initial and final ratio, respectively, 
of the donor star mass to the neutron star mass. Assuming 
the neutron star mass to be constant during the CE-phase 
we can also write k = M2/Mcoie- The value for k is thus 
the mass of the donor star, in units of its core mass, at the 
onset of the RLO. This allows for a convenient rewriting of 
Eq. Q: 



k- 



ao 



1 + 2q{k - l)/77ceAri, 



= [fc + 2go(fc-l)/77ccAri]-i(2) 



The post-CE value for the mass ratio q = Mwd/A/ns — 0.25 
is the present value in PSR J1614— 2230, which is directly 
determined from measurements (see Table[T]). Hence we have 
k — M2/MWD. Taking Mwd ~ 0.50 M0 as the core mass 
we must first determine the value of M2 (and thus k) from 
stellar evolution calculations. To this purpose we used a de- 
tailed one-dimensional hydrodynamic stellar ev olution code. 
This c ode has been described in detail e.g. in iHeger et al.l 
(|2000l ). Using solar chemical abundances {Z = 0.02) and a 
mixing-length parameter of a = l/Hp = 1.5 ('Langeil ll99lh 
we find 2.4 M2/MQ ^ 2.6, see Fig. [T] if we disregard core 
convective overshooting. Including a core conve ctive over- 
shooting parameter of Sqv =0.20 l|Claretl [20071 ) allows for 
donor masses as low as 2.2 Mq to produce a final WD mass 
of 0.50 Mq. Hence, 4.4 < A: ^ 5.2 and we can now use Eq. ((2]) 




Mass ratio parameter, k (q^/q) 

Figure 2. Constraints on stellar parameters assuming common 
envelope evolution in the PSR J1614— 2230 progenitor system. 
The shaded area in the red box indicates the allowed parameter 
space for the radius of the progenitor star at onset of RLO and 
the mass ratio parameter, k. The various curves correspond to 
different A- values of the binding energy of the envelope (see text) . 



to find the pre-CE orbital separation, ao and hence the ra- 
dius of the Roche-lobe filling donor star. 

In Fig. [2] we demonstrate that indeed PSR J1614-2230 
could have evolved from a CE. The shaded rectangular area 
shows the parameter space of solutions. The fc-values are 
constrained by the initial donor mass, M2 which in turn 
is constrained by the observed white dwarf mass, Mwd. 
The upper limit for the radius of the donor star at the 
onset of the RLO is simply its maximum possible radius 
on the AGB, -Rmax. We notice from the curves in the fig- 
ure that only A-values larger than about 2 are in agree- 
ment with this constraint (i.e. of having a donor radius less 
than 7?max). The reason for this is the relatively wide orbit 
of PSR J1614-2230 with an orbital separation of 24.1 ijg. 
Hence, only a modest amount of orbital energy was released 
during spiral-in - almost independent of the pre-CE sepa- 
ration, ao since usually a <^ ao - and therefore the binding 
energy of the donor star envelope cannot have been too large 
for a successful envelope ejection (-Ebind « A~^). The lower 
limit of the progenitor star radius at ^ 300 -R0 is there- 
fore determined by exactly this requirement of having an 
envelope with small enough binding energy (in this case 
corresponding to A ^ 2) such that it can be successfully 
ejected during the spiral-in phase. (For a graphical example 
of a slightly more massive donor star of 3 Mq, see Fig. [1] in 
Dewi & Tauris, 2000). If the donor radius is smaller at the 
time of the onset of the CE then its A-value is too small 
(i.e. its envelope binding energy is too large, on an absolute 
scale, to allow ejection from the available orbital energy re- 
lease). The outcome is a merger event - possibly leading to a 
Thorne-Zytkow object (|Thorne fc Zvtkowiri977h . A similar 
fate is expected for donor stars of late Case B RLO. These 
stars also possess a deep convective envelope, resulting in a 
CE evolution. However, they are less evolved than stars on 
the AGB and have much smaller A-values and hence more 
tightly bound envelopes which strengthens the case for a 
merger. 

To summarize, based on the orbital dynamics and the 
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masses of the two stellar components, Case C RLO (lead- 
ing to a CE and spiral-in) is possible to have occurred in 
PSR J1614-2230. This would have the implication that the 
neutron star was born massive with a mass close to its 
presently observed mass of 1.97 Mq. However, see further 
discussion in Sections [3] and [J] and also in paper II. 



2.1.1 Case BB RLO following a common envelope? 

One might ask whether the 0.50 M© CO WD in 
PSR J1614-2230 could have formed via Case BB RLO from 
a ~ 1 Mq helium star which had previously lost its hydro- 
gen envelope in a common envelope phase. In this case the 
neutron star could have accreted significantly and needed 
not be so massive at birth (only about 1.5 A'Iq if it accreted 
0.5 Mq subsequently). We have tested this hypothesis and 
find that it is not a possible formation channel. The reason 
is that there is not enough orbital energy available (by an 
order of magnitude) to eject the envelope of the WD pro- 
genitor already on the red giant branch (RGB). Assume the 
progenitor of the naked 1 Mq helium core was a 5 — 7 Mq 
star (the exact mass depends e.g. on the assumed amount of 
core convective overshooting). For such a star on the RGB 
we find that the (absolute) binding energy of the envelope is 
1.3 — 2.2 X lO*** erg - according to the formu la for the bind- 
ing en ergy in Section 12.11 and using Fig. 1 in iDewi fc TaurisI 
l|2000h for estimating the largest possible A-values to be 
0.4 — 0.6. The Case BB RLO (from the less massive helium 
star to the more massive neutron star) would have widened 
the orbit to its currently observed orbital period of 8.69 days. 
Assuming a conservative mass transfer the initial orbital sep- 
aration would have been about 10.7 i?0. This value is then 
equivalent to the post-CE orbital separation corresponding 
to an orbital energy of —2.7 x 10*^ erg. Hence, even if all 
available orbital energy was released from an initial separa- 
tion at infinity there would be far too little energy to eject 
the envelope of the RGB star (Ai5orb "C -Ecnv). Only on the 
AGB is the envelope loose enough to allow ejection from 
spiral-in of the neutron star. But in this case the core has 
already burned its helium into carbon and oxygen and then 
we are back where we began in Section [2TT] 



2.2 Early Case B RLO - thermal timescale mass 
transfer 

If a 3 — 5 Mq donor star fills its Roche-lobe shortly af- 
ter leaving the main sequence {early case B) its enve- 
lope is still radiative and the binary may survive ther- 
mal timescale mass transfer. T his was shown a de cade 
ago in three independent papers: King &: RitteJ (|l999f ) and 
iPodsiadlowski fc RappaportI (|2000') studied the formation 
and evolution of Cyg X-2 and Tauris ct al. (2000) inves- 
tigated the formation of binary millisec ond pulsars with 
a CO WD compani on. Ahhough .Tauris et al.l l|2000l ) and 
IPodsiadlowski et al.l (|2002l ') have demonstrated that one can 
form systems with a 0.50 Mq CO WD and Poib = 8.7 days 
(as observed in PSR J1614— 2230) they both assumed in 
their calculations an initial canonical neutron star mass of 
1.30 — 1.40 Mq, which does not apply in this scenario since 
the neutron star only accretes a few 0.01 Mq during the 
Case B rapid mass-transfer phase, thus disqualifying the 



neutron star from reaching its present mass of 1.97 Mq. If 
the initial mass ratio between the donor star and the neutron 
star is of the order go — 2 — 3 the orbit shrinks significantly 
in response to mass loss, as mentioned previously. This leads 
to highly super-Eddington mass-transfer rates and hence we 
can apply the isotropic re-emission m ode of mass transfer 
iBhattacharva fc van deri Heuvej|l99lh . In this model mat- 
ter flows over from the donor star (M2) to the accreting neu- 
tron star (A/ns) in a conservative manner and thereafter a 
certain fraction, /? of this matter is ejected from the vicinity 
of the neutron star with the specific orbital angular momen- 
tum of the neu tron star (for example, in a jet as observed in 
SS433, see also lKing fc Begelmanlll999l ) . Integrating the or- 
bital angular momentum balance equation one can find the 
change in orbital separation during the is otropic re-emission 
RLO fe.g. lTauris|[l99^ : lKing et al.ll200ll ): 



a 

0.0 



go(i-/?) + i 
g(i -/?) + ! 



3/3-5 
1-/3 



go + 1 



(3) 



where it is assumed that P remains constant during the 
mass-transfer phase. Indeed iTauris et al.l (|2000l ) showed in 
their Fig. 1 that in intermediate-mass X-ray binary (IMXB) 
systems very little mass is accreted onto the neutron star 
since the timescale for the mass-transfer phase is very short 
(~ 1 Myr) leading to highly super-Eddington mass-transfer 
rates by 3-4 orders of magnitude and hence (3 > 0.999. It is 
therefore interesting to con sider Eg. ll3l) in the limit where 



/3 — >■ 1 and we find (see also lKing et al. 
in orbital period: 



lim 



till] fc^e^'fi-*' 

g + 1 



for the change 



(4) 



under the above mentioned assumptions and by applying 
Kepler's third law. 

In Fig.|3]we demonstrate that early Case B mass trans- 
fer is not possible to have occurred in the progenitor binary 
of PSR J1614— 2230. The constraints on k for Case B mass 
transfer can be found from Fig[T] For Case B mass trans- 
fer (between evolutionary epochs TAMS and the RGB) a 
progenitor star of 4.0 — 4.5 Mq is needed to yield a core 
mass of 0.50 AIq (the observed mass of the white dwarf in 
PSR J1614-2230). Hence, we find 8 < A; < 9 for this sce- 
nario. Recalling q — Mwd/Mns = 0.25 and given the or- 
bital period of P = 8.69 days then, according to Eq. (Q, 
this would require an initial orbital period, Po — 0.7 days 
which is not possible for Case B mass transfer - the mini- 
mum initial period for (early) Case B RLO is shown as the 
red line in Fig. |31 In fact with such a short initial orbital pe- 
riod the donor star would fill its Roche-lobe radius instantly 
on the ZAMS. Even if we expand the interval of donor star 
masses to the entire range 2.5 < M2/MQ < 6.0 early Case B 
RLO would still not be possible to explain the parameters of 
PSR J1614— 2230. We can therefore safely rule out all vari- 
ations of Case B mass transfer (early, late and Case BB). 



2.3 Case A RLO — mass transfer from a main 
sequence star 

In order to reproduce PSR J1614-2230 via Case A RLO we 
notice from Fig[T]that we must at first glance require an ini- 
tial donor mass of almost M2 — 5 Mq in order to end with 
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Moss ratio parameter, fc (q^/q) 

Figure 3. The blue line represents the initial orbital period of the 
progenitor X-ray binary as a function of the mass ratio parameter 
k using the isotropic re-emission model for early Case B RLO. The 
two vertical dotted lines indicate the interval of possible values of 
k for PSR J1614-2230. The hatched region is excluded for Case B 
RLO since the donor star would have filled its Roche-lobe before 
reaching the stage of shell hydrogen burning. The original donor 
star mass is M2 = k X 0.50 M©. Note, the fe- values for Case B 
RLO are larger than the fc-values expected for Case C RLO in 
Fig. [2] The reason for this is that the core of a Case B donor 
star has not had time to evolve to the large core masses found 
in Case C RLO (see Fig. [TJ , resulting in higher required donor 
masses, M2 for Case B compared to a Case C scenario. 



a final white dwarf mass of 0.50 M©. However, tiie evolu- 
tion of Case A RLO is somewhat compl ex and not straight 
forward to analyse analytically (see Tauris fc LangeJ[201l] . 
for further details). The estimated TAMS core masses from 
Fig[l]are not necessarily good indicators for the final mass 
of the white dwarf remnants evolving from Case A donors in 
X-ray binaries for two reasons: 1) forced mass loss from the 
Roche-lobe filling donor star results in a lower core mass 
as the donor now evolves less massive, and 2) the forma- 
tion of an outgoing hydrogen shell source during the final 
phase (phase AB, see below) of the mass transfer causes 
the core mass to grow with the helium ashes left behind. 
Therefore, to obtain the final mass of the white dwarf re- 
quires detailed numerical stellar models. The overall effect 
is that the core mass will have grown somewhat by the time 
the system detaches from the RLO. Hence, the white dwarf 
remnant left behind is expected to be slightly more massive 
than the donor core mass at the TAMS. For this reason the 
ZAMS mass interval found from Fig. [1] for a Case A donor 
star of PSR J1614— 2230 should be considered as an upper 
limit and in the following we explore donor masses down 
to 4.0 Mq. (Below this limit the WD remnant becomes too 
Ught, see Section l2.3.ip . Stars more massive than 5.0 
could leave behind a core mass even less than 0.50 Mq if the 
mass transfer is initiated well before reaching the TAMS. 
However, these binaries would not be dynamically stable 
with a neutron star accretor, see Section 

Our analysis reveals the parameter space of Case A 
binaries which produce the characteristic parameters of 
PSR J1614— 2230. Figs. 4-8 show an example of a calcu- 



lation of a possible progenitor X-ray binary. This IMXB 
started out with an initial donor star of mass 4.50 Mq and 
a 1.68 A/q neutron star accretor having an initial orbital 
period of 2.20 days. We assumed here a convective core- 
overshooting of 5ov — 0.20. In Fig. |4] and Fig. [5] we demon- 
strate that the system experiences three phases of mass 
transfer (hereafter denoted phases Al, A2 and AB, respec- 
tively) . In phase Al the mass transfer proceeds on the ther- 
mal timescale (see lLanger c t al. 2 000i , for further discussions 
on thermally unstable mass transfer). The reason for this 
is the initially large mass ratio between the heavier donor 
star and the lighter neutron star which causes the orbit to 
shrink in response to mass transfer. As mentioned earlier, 
the outcome in this situation is that the donor star over- 
fills its Roche-lobe even more - leading to further mass loss 
- and within 1 Myr it looses more than 3Mq at a rate 
exceeding 10~^ Mq yr^^. This rate is still less than the esti- 
mated limit for which photons are trapped, leading possibly 
to rapid neutrino cooling and hyper critical accretion. I t has 
been demonstrated, for example bv lFrver et all (|l996l ) and 
iKing fc BegelmanI (Il999l ). that it takes an accretion rate of 
at least a few times lO~'*M0yr~^ before the outward dif- 
fusion of photons cannot overcome the inward advection of 
photons in the accreted matter. (The exact limit is uncertain 
and depends, for example, on the amount and the geometry 
of outflows). Although the donor star is driven out of ther- 
mal equilibrium during phase Al it manages to retain hydro- 
static equilibrium and the system can in this case avoid a so- 
called delayed dynamical instabi lity (iHiellming fc WebbinkI 
ll987l : lKalogera fc Webbinkll996h which would have resulted 
in a common envelope and most likely a merger event. 

The final mass of the neutron star in our example is 
1.99 Mq. The neutron star has thus accreted a total of 
0.31 Mq. The amount accreted in each phase is found from 
Fig. |4] by integrating the area under the blue line which 
falls below the Eddington accretion limit (red dashed line). 
Hardly any accretion takes place during the very short (ther- 
mal timescale), ultra super-Eddington phase Al. Phases A2 
and AB proceed on nuclear timescales dictated by core burn- 
ing of the remaining hydrogen and, later on, hydrogen shell 
burning, respectively. 

Fig. [6] shows the track of this IMXB donor star in 
the HR-diagram on its path to forming a carbon-oxygen 
white dwarf (CO WD) orbiting a millisecond pulsar. The 
Case A RLO mass transfer is initiated at an orbital period of 
2.20 days. At this stage the core of the donor star still has a 
central hydrogen mass abundance of Xc = 0.09. The error of 
putting the donor star on the ZAMS is small given that the 
progenitor star of the neutron star (i.e. the primary star of 
the ZAMS binary) most likely had a mass of at least 20 Mq 
(see arguments in Sections |3] and |4ll and hence a lifetime 
of less than 10 Myr, which is short compared to the main 
sequence lifetime of a 4.50 Mq star. Binaries with shorter 
initial periods will have less evolved donor stars when en- 
tering the mass-exchange phase. This leads to lower helium 
cores masses which are then often below the threshold for 
igniting the triple-a process. Hence, these systems will leave 
behind pulsars wit h a low-mass helium WD c ompanion, as 
first pointed out bv lPodsiadlowski et all tOOj) - for further 
discussion on these systems, see Tauris fc Langej (|201ll ). 

When the donor star settles on the He-ZAMS its lumi- 
nosity is entirely generated by the triple-a process in the 
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Figure 4. Case A RLO in a binary initially made of a 4.50 Mq 
donor star with a 1.68 Mq neutron star and an initial orbital 
period of Porb = 2.20 days. The graph shows the mass-transfer 
rate from the donor star as a function of its age. Three phases of 
mass transfer (Al, A2 and AB - see text) are identified. 









1 ' 1 ' 1 ' 1 
4.50 Mq donor star 


OV=0.20 


1 ' 1 ' 




1 .68 Mq neutron star 








Pj^jt = 2.20 days 


















RLO Case A 








(X^=0.09) 






% \ 

w 


1 





5 4.8 4.6 4.4 4.2 4 3.8 



i°g (K) 

Figure 6. Evolution of the mass loosing donor star in the HR- 
diagram. Starting and termination points of the three phases of 
mass transfer are shown by filled circles, open triangles and open 
squares, corresponding to phases Al, A2 and AB, respectively (cf. 
Fig- Ell - The core helium burning is ignited after detachment from 
the Roche-lobe, see text. 
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Figure 5. The mass-transfer rate as a function of the decreasing 
donor star mass for the stellar evolution calculation shown in 
Fig. |4] Very little mass (~0.01Mq) is accreted by the neutron 
star during phase Al which proceeds on a thermal timescale. 



core - marked in Fig. |6] with the symbol "3a". At this 
stage the donor star has an age of 180 Myr. The ignition 
of the triple-Q process actually occurs gradually already 
from (logTeff, log(L/LQ)) = (4.2, 2.9) when Tcorc > 10* K, 
shortly after the detachment from the Roche-lobe, and over- 
laps with the ceasing stages of the shell hydrogen burning. 
The curly loop at \og{L / Lq) ~ 2 indicates the beginning of 
the shell helium burning phase at an age t — 300 Myr. 

The orbital evolution is shown in Fig. [7] where the or- 
bital period is plotted as a function of decreasing donor star 
mass. The final orbital period of our system is 8.67 days. The 
widening of the orbit is quite significant in phase AB where 
the mass ratio, q is small. It is also during phase AB that the 
neutron star (NS) gains the majority of its accreted mass - 
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Figure 7. Orbital period (red line) and pulsar mass (green line) 
as a function of decreasing donor star mass. The symbols are 
equivalent to those defined in Fig. |6l 



see green line in Fig. [T] The donor star of the X-ray binary 
Cyg X-2 is an example of a hy drogen s hell burning star nea r 
the end of phase AB l|Podsi adlowski fc RappaportI |2000| ) . 
Interestingly enough, a massive neutron star (~ 1.8 M©) 
has been infe r red fo r this source as well, see for example 
ICasares et al] (ll998D . 

The chemical abundance profile of the resulting CO WD 
is shown in Fig.|8] We notice that the inner core (~ 0.28 Mq) 
contains almost 90% oxygen (mass fraction). The CO WD 
is seen to have a 0.04 Mq helium envelope and a tiny 
content of up to 8% hydrogen in the outermost lO"'* Mq 
(amounting to a total of 1.7 x 10~^ Mq). The oxyg e n con - 
tent of the WD is larger than found by iLin et al] (1201 ll ) 
using the MESA code (see their Fig. 6, panel d). To further 
test the chemical composition of our WD we compared a run 
of a similar binary system (4.50 Mq donor star, 1.68 Mq NS, 
Porb ~ 2.20 days, Sqv — 0.20) using the Eggleton stellar 
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Figure 8. The chemical abundance structure of the CO WD 
formed in the Case A RLO shown in Figs. I4I7I This profile 
could well resemble the structure of the CO WD companion to 
PSR J1614— 2230 — possibly applicable to a WD cooling model. 
The radius of the WD, once it settles on the cooling track, is 
9500 km yielding a surface gravity of log(g) = 7.9. 



evolution code. This computation yielded an oxygen con- 
tent in agreement with our result (Evert Glebbeek, private 
co mmunicat i on). T he reason for the different oxygen yield 
m iLin et al.l (|201ll ) could very well be caused by the use of 
different ^^C(a.,^Y^O reaction rates. 

To summarize, the final outcome of our example shown 
for Case A evolution is a 0.501 A/q CO WD orbiting a 
1.99 Mq (millisecond) pulsar with an orbital period of 
8.67 days - almost exactly in agreement with the observed 
parameters of PSR J1614-2230, see Table [T] 

2.3.1 Permitted parameter space for Case A RLO leading 
to the formation of PSR J 1614- 2230 

We have demonstrated above that both Case C and Case A 
RLO during the X-ray phase can reproduce the observed 
parameters of PSR J1614— 2230 for suitable initial masses 
of the two components and their orbital period. In order 
to search the entire parameter space of Case A systems we 
explored a range of binaries by altering the stellar masses, 
the orbital period and the accretion efficiency. In Fig. [9] we 
show the grid of resulting NS+WD systems in the final or- 
bital period versus final neutron star mass plane. This plot 
was obtained by varying the initial mass of the neutron star 
as well as the accretion efficiency for a fixed value of the 
donor star mass, M2 ~ 4.50 Mq at the onset of the X-ray 
phase. It is obvious that the final mass of the neutron star is 
a growing function of its initial mass as well as the efficiency 
of accretion. 

In order to be able to compare with previous work we 
have in this plot defined the accretion efficiency as a value 
in percent of the Eddington mass accretion limit (AfEdd) for 
pure hydrogen on a neutron star with a radius of 10 km, such 
that a value of 100% corresponds to the canonical accretion 
rate of 1.5 x 10"® Mq yr"^ A value larger than 100% cor- 
responds to either accretion at slight super-Eddington rates 
or accretion of matter with a larger mean molecular weight 



per electron (e.g. an accretion efficiency value of 200% cor- 
responds to accreting pure helium at the Eddington limit). 
Many of the grid points in Fig. are not obtained from 
actual stellar evolution calculations. For example, the evo- 
lution leading to grid points based on an initial neutron 
star mass of 1.4 Mq were dynamically unstable in our mod- 
els, leading to runaway m ass transfer (see below, and also 
IPodsiadlowski et al.! '2002'). Neverthele ss, one can st i ll com - 
pare with the calculations in .Podsiadlowski et al. I (|2002l ). 
Using a neutron star with an initial mass of 1.4 Mq orbiting 
a donor star of mass 4.5 Mq with an initial orbital period of 
2.38 days, and assuming an accretion efficiency of 50%, these 
authors end up with a NS-l-WD binary with a final neutron 
star mass of 1.507 Af©, a white dwarf mass of 0.471 Mq and 
an orbital period of 3.43 days. (Also in their work they find 
that such an X-ray binary is barely on the edge of stability 
and note that this system may be dynamically unstable). 
The result of their calculation is shown in our figure with 
an open black circle and the agreement with our result is 
indeed quite good (cf. the orange neighbour point in our 
grid just below their point). Our result is based on one of 
our calculations with a 4.50 Mq donor star and a neutron 
star mass high enough to avoid a dynamical instability, for 
example of mass 1.7 Mq. The effect of changing the neutron 
star mass and/or the accretion efficiency can easily be found 
analytically from an extrapolation of our calculated model 
using Eq. ^ for each of the three phases of mass transfer 
by adapting the new values of /3, g and qo. The underly- 
ing assumption that the amount of mass transfered from 
the donor star remains roughly constant (i.e. independent 
on the neutron star mass) has been tested by us and shown 
to be correct. This was done by directly comparing the re- 
sult of an extrapolated model with a calculated model. See 
iTauris fc Langeij l|201ll) for further discussion. 
In our stellar evolution code the Eddington accretion limit 
(i.e. the accretion efficiency) depends on the chemical com- 
position of the accreted matter as well as the radius of the 
neutron star, both of which are time dependent - see de- 
scription in Paper II, and the red dashed lines in Figs. 4-5. 

It is important to notice from Fig.|9]how the final orbital 
period is correlated with the initial mass of the neutron star 
(increasing in value upwards in the grid diagram, see solid 
lines). The grid clearly shows that, in the frame of Case A 
RLO, the neutron star in PSR J1614— 2230 cannot have been 
born with the canonical birt h mass of about 1.3 Mq. This 
conclusion was also found bv iLin et al.l (|201ll ). 

Using our stellar evolution code we find that our IMXBs 
are only stable against dynamical mass transfer for initial 
mass ratios up to go — 2.7 — 3, e.g. corresponding to initial 
donor masses at most 3.5 — 4.0 Mq for a 1.3 Mq neutron star 
and (what is important for PSR J1614— 2230) donor stars 
up to 5.0 AfQ for a 1.7 A^q neutron star. Therefore we adapt 
5.0 Mq as the upper limit for the initial mass of the donor 
star. The lower limit for the mass of the donor star is con- 
strained by the mass of the CO WD in PSR J1614-2230. We 
find a lower limit of about 4.0 Mq (a 3.5 Mq donor star in 
the region of relevant initial orbital periods leaves behind a 
WD mass of only 0.39 A/q which is 20% smaller than needed 
for PSR J1614-2230). 

The effect on the final orbital period and neutron star mass, 
imposed by changing only the donor star mass and keep- 
ing all other parameters fixed, can be visualized by mov- 



PSR J1614—2230: Evidence for a NS horn massive 9 



~o to 



- 


1 ' 1 ' 1 ' 

^= 

M2 = 4.50 Mg "o 


o 

Q 


- 




1.8 Mg 


,o 






1 .7 Mg . . 


J 














1.6 Mg . 














I.I.I, 





1.2 1.4 1.6 1.8 2 2.2 2.4 



Final neutron star mass (M/Mg) 

Figure 9. The final orbital period as a function of final neutron 
star mass for a grid of X-ray binaries evolving from a 4.50 Mq 
donor star through Case A RLO. In all models the CO WD is 
formed with a mass of about 0.51 ± 0.01 M0. The initial orbital 
period was in all cases about 2.2 days, corresponding to a core 
hydrogen content of ~ 10% at the time of RLO. The variables 
are the initial neutron star mass (solid lines) and the accretion 
efficiency (dotted lines). The observed values of PSR J1614— 2230 
are shown with a red star. Our calculations show that indeed 
PSR J1614— 2230 could have evolved from a 4.5 Mq donor star 
and a neutron star born with a mass of ~ 1 . 7 Mq , accreting at 
an efficiency of 150% — see text. 



ing the entire grid in Fig. [§] up or down for a less massive 
and a more massive donor star, respectively. We find that 
a 4.0 Mq donor star would need to be in a binary with a 
neutron star of initial mass of 1.55 M© in order to repro- 
duce PSR J1614-2230. The 5.0 Mq donor star would need 
a neutron star of initial mass of 1.77 Mq in order to repro- 
duce PSR J1614— 2230. In both cases the required accretion 
efficiency value is about 160-170%. However, the precise lim- 
its depend on, for example, the uncertain strength and the 
assumed underlying physics of the tidal torques and result- 
ing spin-orbit coupUngs (which may help to stabilize the 
orbital evolution in X-ray bin aries with even higher mass 
ratios, iTauris fc Savoniiell200il ) . 

There is evidence from comparison of numerical binary 
stellar evolution calculations and observations of recycled 
pulsars that the accretion efficien cy in som e cases is rather 
low (see Fig. 4b, §5.7 and §6.4 in iTauris fc Savoniic 1999), 
possibly as a result of the propeller effect and/or accretion 
disk instabilities. A low accretion efficiency would have the 
implication for PSR J1614— 2230 that it was born with an 
even higher mass. For a 4.5 Mq donor star and assuming 
an accretion efficiency of 50% (according to the definition 
given earlier) this would result in a neutron star birth mass 
of about 1.9 Mq. However, in this case one can only repro- 
duce the observed orbital period of 8.69 days by increasing 
the specific orbital angular momentum of the lost material 
beyond the value expected for matter in the vicinity of the 
accreting neutron star (i.e. by matter ejected from a loca- 
tion elsewhere within the binary system). To illustrate this 
one can see from Fig. [9] that an initial neutron star mass of 
1.9 Mq and an accretion efficiency of 50% would yield a final 



orbital period of about 18 days (much larger than observed 
in PSR J1614-2230). This problem can be solved if one 
instead assumes that the material lost from the system has 
the average specific orbital angular momentum of the binary. 
The reason is that the orbit widens mainly during the mass- 
transfer phase AB, when the accreting neutron star is more 
massive than the donor star, and in our model the ejected 
material has the low specific orbital angular momentum of 
the neutron star during this phase. 



3 EVOLUTION OF PROGENITOR BINARIES 
FROM THE ZAMS TO THE X-RAY PHASE 

In the previous section we presented evidence for two dif- 
ferent formation scenarios (hereafter simply called Case A 
and Case C) for the formation of PSR J1614-2230. Hence, 
we know the required parameters at the onset of the RLO, 
for the X-ray binary containing a neutron star and a non- 
degenerate star, and one can then try to calculate backwards 
and estimate the initial configuration of ZAMS binaries 
which may eventually form a system like PSR J1614— 2230. 
Our brief description presented here is only qualitative. An 
analysis including, for example, dynamical effects of asym- 
metric supernovae (SNe) is rather cumbersome. To obtain 
a set of more detailed and finetuned parameters of the pro- 
genitor binaries one would need to perform a population 
synthesis (which is beyond the scope of this paper). Nev- 
ertheless, below we present the main ideas. The results of 
our simple analysis are shown in Table [2] and illustrated in 
Fig.m 

A mass reversal between the two interacting stars dur- 
ing the evolution from the ZAMS to the X-ray phase can 
be ruled out in both scenarios as a consequence of the large 
difference in ZAMS mass between the two stars. The sec- 
ondary star, the donor star of the X-ray binary, had a mass 
of 2.2 — 2.6 M0 or 4.50 Mq, respectively. The primary star 
had at least a mass above the threshold of ~ 10 Mq for pro- 
ducing a NS in a close binary. In Section [4.21 we argue that 
the mass of the progenitor of the NS in PSR J1614-2230 
was most likely more than 20 Mq since it left behind a mas- 
sive NS. Hence, there is no doubt that the NS is the rem- 
nant of the original primary star. Mi (i.e. the initially more 
massive) in the ZAMS binary and the donor star in the 
X-ray phase thus descends from the secondary ZAMS star, 
M2. Based on this argument, due to the small mass ratio 
M2/M1 ~ 0.1 — 0.25 we would expect the progenitor binary 
to have evolved through a CE and spiral-in phase in both 
Case A and Case C on their path from the ZAMS to the SN 
stage. The initial ZAMS orbital period was then probably 
quite large (> 10^ days) in order to let the primary star 
evolve to a late evolutionary stage before the onset of the 
CE. This would both facilitate the ejection of the envelope 
- which is more loosely bound at late evolutionary stages - 
and ensure that the helium core of the neutron star progen- 
itor evolved "clothed" - see Section r4.2l We will now briefiy 
discuss each of the two scenarios: Case C and Case A. 

3.1 Case C progenitor binary 

The aim here is to obtain a 2.2 — 2.6 Mq non-degenerate 
AGB star orbiting a NS (see Section [2. Hence, the binary 
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Table 2. Evolution and characteristics of the possible progenitor binaries of PSR J1614— 2230. The different columns correspond to the 
different cases of RLO in the X-ray binary phase. In this table the evolution with time goes from top to bottom (see also Fig. I10|l . 
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Figure 10. An illustration of the progenitor evolution leading to 
the formation of PSR J1614-2230 for both Case A and Case C. 
Only a few evolutionary epochs are shown for simplicity. 



must have been very wide (~ 10'^ days) following the SN in 
order to allow the donor star to ascend the AGB before ini- 
tiating mass transfer. A wide post-SN orbit could have been 
the re sult of a large kick im parted to the newborn NS at 
birth (|Lvne fc Lorimeilll994 ). However, tidal interactions in 
a highly ec centric orbit ma y significantly reduce the orbital 
separation (|Sutantvdll97i ;) thus preventing the companion 
star from ascending the AGB before it fills its Roche-lobe. 



The other alternative, namely that the wide post-SN orbit 
simply reflects that the pre-SN orbit was wide too is also pos- 
sible, but not very likely either. Such a wide system would 
rarely survive any kick imparted to the newborn neutron 
star - and a small kick originating from an electron capture 
supernova would be in contradiction with the hi gh neutron 
star mass (Podsiadlowski et al ] |2004l :l van den Heuvel 2004). 
Furthermore, the expected outcome of the first CE phase 
(between the ZAMS and SN stages) is not a wide system. 
The reason is that the binding energy of the massive primary 
star's envelope is too large, on an absolute scale, to allow for 
an early ejection. Hence, there could probably not have been 
a "mild in-spiral" as a result of an easy ejection of the enve- 
lope resulting in only a modest conversion of orbital energy 
and allowing the orbit to remain fairly wide. An alterna- 
tive possibility is that the initial orbit was so wide that the 
two stars did not exchange mass during the giant phase of 
the primary star. The subsequent kick in the SN then shot 
the newborn NS into a closer orbit around the secondary 
star (JCalogera 1998). However, this scenario would require 
a very fortunate finetuning of both the kick magnitude and 
the direction, making it unlikely too. 



3.2 Case A progenitor binary 

To produce an X-ray binary with an orbital period of only 
2.20 days and a 4.50 Mq donor star seems more likely com- 
pared to the case described above. The short orbital period, 
both before and after the SN, is a simple consequence of the 
in-spiral during the first CE-phase when the primary star 
was a giant. We therefore conclude at this stage, that based 
on binary evolution considerations PSR J1614— 2230 seems 
more likely to have evolved from a Case A RLO X-ray bi- 
nary and that the initial ZAMS system was composed of a 
^ 20 Mq primary with a 4 — 5 Mq secondary in a wide orbit. 
In the next section we discuss the stellar evolution of the NS 
progenitor in much more detail and in Paper H we discuss 
progenitor systems based on the spin-up of the neutron star. 
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4 DISCUSSION 

The precise measurement of the high neutron star mass 
in PSR J1614— 2230 leads to interesting implications for 
the nuclear physics beh ind the equation-of-state (e.g. 
iLattimer fc PrakashI [ioiol '). Equally important, the result 
has caused renewed interest in modelling close binary evo- 
lution - in particula r the m ass-transfer phase. In a recent 
paper iKiziltan et al.l (|201lh discuss the possibility of " al- 
ternative evolution" in order to explain the observed mass 
of PSR J1614-2230. However, in Section [2] of this pa- 
per we have demonstrated that this is not required and 
PSR J1614— 2230 may have followed standard evolution 
paths expected from stellar astrophysics. 
The challenge in reproducing massive binary millisec- 
ond pulsars with a CO white dwarf companion (like 
PSR J1614— 2230) is to get all three fundamental observable 
parameters correct: the masses of the two compact objects 
and the orbital period. In this paper we have demonstrated 
a methodical approach to do this involving all three mass- 
transfer scenarios (RL O Cases A, B an d C). 

In a recent paper iLin et all (|201ll ) systematically com- 
puted the evolution of a large number of Case A and early 
Case B IMXB binaries and applied their results to under- 
stand the formation of PSR J1614-2230. They conclude 
that a system like PSR J1614— 2230 requires a minimum 
initial neutron star mass of at least 1.6 ± 0.1 Mq, as well as 
an initial donor mass of 4.25 ±0.10 Mq and orbital period of 
^49 ± 2 hr (2.05 ±0.1 days). In general their Case A results 
are in fine agreement with our Case A results. The main dif- 
ference is their rather narrow range of required donor star 
masses (4.25 ± 0.10 A/0) compared to our wider interval of 
4.0-5.0 Mq. This minor discrepancy could arise from using 
different values of the convective core-overshooting parame- 
ter, the mixing length parameter and/or the chemical com- 
position of the donor stafl However, it is interesting to no- 
tice the broad agreement in the final results given that the 
stellar evolution codes are different. 

It is important to emphasize that other observed re- 
cycled pulsars with a CO WD companion and Porb — 
5-15 days (such as PSR J0621±1002 and PSR J2145-0750) 
are, in general, not expected to have a massive neutron star. 
iTauris et all ( |200d !) demonstrated that these systems can be 
reproduced b y early Case B RLO wi th a typical 1.3 Mq neu- 
tron star and lvan den Heuvell (| 19941 ) argued for a CE evolu- 
tion scenario using a neutron star mass of 1.4 M©. Systems 
with a massive neutron star, like PSR J1614— 2230 inves- 
tigated here, require the neutron star to be born massive 
- possibly followed by an extended phase of mass transfer 
allowing for significant further accretion of matter. 



4.1 Neutron star birth masses in pulsar binaries 

The interval of known radio pulsar masses ranges from 
1.17 Mq in the doub le neutron star binary PSR J1518±4909 
(3-(T upper limit, IJanssen et al.l |2008| ) to 1.97 A/0 in 
PSR J1614— 2230, discussed in this paper. The most massive 
of the non-recycled companions in double neutron star sys- 
tems is the unseen companion in PSR 1913±16 which has a 



^ These parameters are not stated in their present publication. 



mass of 1.389 A/0 JWeisberg et al.ll2010l V Interestingly, the 
observed pulsar in this binary is the most massive of the 
(mildly) recycled pulsars detected in any of the ten dou- 
ble neutron star systems. It has a mass of 1.440 A/q. How- 
ever, the relatively slow spin period of this pulsar (59 ms) 
hints that only about 10~^ Mq was needed in the recycling 
process (see Paper II) and thus 1.44 Mq is the previously 
known upper limit derived for the birth mass of any neu- 
tron star detected in a binary pulsar system. Only a few of 
the ~ 120 binary pulsars with WD companions have mea- 
sured masses - see Paper II - and just a handful of these 
are more massive than 1.44 Mq. But even in those cases the 
mass determinations are often very inaccurate and also in- 
clude the mass accreted from the progenitor of their WD 
companion. However, in this paper we have demonstrated 
that the birth mass of the neutron star in PSR J1614— 2230 
is at least 1.7±0.15 A/0. This result is important for under- 
standing the physics of core collapse supernova and neutron 
star formation. 



4.2 Predictions from stellar evolution and 
SN explosion models 

The difficulty of predicting an upper limit for the birth 
mass of a neutron star is mainly caused by unknown de- 
tails of stellar evolution and explosion physics, as well as the 
neutron star equation-of-state. For example, internal mixing 
processes, wind mass loss and key nuclear reaction rates are 
still not known accurately. Additionally, even models with 
the same input physics predict a rather erratic behaviour of 
the final iron c ore mass as a function of ZA MS mass - cf. 
Fig. 1 7 in lWooslev et all (|20o3 ') and Fig. 4 in lTimmes et all 
(jl996t ). Moreover, rotation and metallicity (and perhaps 
magnetic fields) affect the final remnant mass too. Rotation 
induced chemical homogeneous evolution will ca use rapidly 
rotating stars to develop larger cores (e.g. lYoon et al. 200a ) 
and a low metallicity content also causes stars to develop 
some what larger core s as a resu lt of reduced wind mass 
loss JVink et a .IbOOlT) - see also iLinden et aD (|2010l ) and 
iMirabel et all l|201ll ) who discuss the relatively high num- 
ber of black hole HMXBs at low metallicities (compared to 
HMXBs with a neutron star companion). 

As for the explosion physics there are, for example, un- 
certainties in deriving the explosion energy as function of 
the pre-supernova structure, and in the determination of 
the mass cut separating the initial compact remnant from 
the ejected matter. Also the amount of fall back material is 
uncertain. 

Despite the above-mentioned uncertainties we can iden- 
tify two main factors which determine the remnant mass of 
a given early-type star: 1) its ZAMS mass, and 2) whether 
or not it looses its hydrogen-rich envelope (e.g. as a result 
of mass transfer in a binary) before or during core helium 
burning. Both of these factors infiuence the carbon abun- 
dance at core helium exhaustion which plays a crucial role 
for the subsequent carbon, neon and oxygen burning stages 
which again determines the size of the silicon and iron cores 
(and thus the mass of the newborn neutron star). 
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4-2.1 To burn carbon convectively or not - the role of the 
^^C/^^O-ratio at central helium depletion 

Following iBrown et all lj200ll ). the ^^C/^^'O-ratio at cen- 
tral helium depletion is determined by the competition 
between the formation of carbon via the triple-Q process 
(3a —> '^^C) and the destruction of carbon, mainly via a- 
capture: ^■^C(q, 7)^^0. More massive stars perform core he- 
lium burning at higher temperatures and lower densities, 
which decrease the net carbon yield. In stars initially more 
massive than about 19 — 20 Mq the resulting core carbon 
abundance is too small (in mass fraction Cc ^ 0.15) to pro- 
vide a long lasting conv ective core carbon burning phase 
(jWooslev fc Weaver! I995I ). The short lasting (radiative) core 
carbon burning phase results in less energy carried away by 
neutrinos. As a consequence, the entropy in the core remains 
fairly high which leads to more massive pre-SN cores. A low 
carbon abundance also causes the carbon burning shells to 
be located further out which increases the size of the carbon- 
free core and thus increases the mass of the iron core to 
form. Single stars with ZAMS masses less than 20 AIq , on 
the other hand, deplete core helium burning with a relatively 
high ^^C/^^O-ratio. These stars undergo significant convec- 
tive carbon burning which subsequently leads to relatively 
small iron cores and hence low-mass neutron stars. 



4-2.2 Single stars / wide orbit binaries (Case C RLO) 

Based on the core c a rbon burning dichotomy discussed 
above, iTimmes et al.l (|l996l ) found a bimodial distribu- 
tion of neutron star birth masses with narrow peaks at 
1.28 ± 0.06 M© and 1.73 ± 0.08 M© (gravitational masses), 
for single star progenitors (type II SNe) below and above the 
critical ZAMS mass of ~ 20 Mq . Based on assumptions of 
rathe r soft equations-of-s tate for nuclear matter at high den- 
sities |Brow^|er|al] (|200lh did not predict high- mass neutron 
stars (~1.7M0) to be formed and therefore concluded that 
all single ZAMS stars in the interva l 20 — 25Mq end their 
lives as black holes - see also iFrveJ (| 19991 . '200^V However, 
given the discovery of the massive pulsar PSR J1614— 2230 
this picture has to change. Following our work in this paper 
the birth mass of this neutron star is at least 1.7 ± 0.15 Mq 
and given the expected wide orbit of its progenitor star (see 
Section[3]) we conclude that it formed from an original ZAMS 
star of mass 20 — 25 Mq which underwent Case C RLO (lead- 
ing to the CE- phase prior to the SN - see Fig. llOp . Hence, 
in general we expect 20 — 25 Mq single (or wide orbit stars) 
to form high- mass neutron stars. Beyond 25 M©, the large 
binding energies of the mantle is expected to result in large 
fall back and the production of a black hole. 

For the sake of completeness, at the other end of the 
scale where stars with an initial mass of about 10 Mq end 
their lives the minimum gravitational mass expected for a 
newborn neu tron star is about 1.25 Mq for an electron cap- 
ture SN fe.g. lPodsiadlowski et al.ir2004l 'l. However, for single 
stars, the initial mass range for stars to undergo electron 
capture SN is very small, such that only a few percent of 
all core c ollapse supernova e are expected to go through this 
channel (|Poelarends et aL| [2008). It is interesting to notice 
that small iron core collapse SNe of type II als o seem to allow 
for th e formation of ~ 1.15 Mq neutron stars (|Timmes et al.l 
1 19961 ). 



To summarize, one would expect three peaks in the dis- 
tribution of initial neutron star masses from the evolution 
of single stars or stars in wide orbit binaries (as pointed out 
bv Ivan den Heuvell 120041 '): one (smaU) peak at ~ 1.25 
from electron capture SNe of 8 — 10 Mq stars (leading to a 
small kick), one peak at 1.25 — 1.4 from iron core col- 
lapse SNe of 10 — 20 Mq stars and one high-mass peak at 
^ 1.7 Af© from iron core collapse SNe of 20 — 25 stars. 
Above a ZAMS mass of 25 Mq single stars are expected to 
form black holes - unless the metallicity is very high (above 
solar). Rapid rotation, on the other hand, may produce black 
holes from stars less massive than 25 Mq- 

4-2-3 Close binary stars (Case A and Case B RLO) 

The progenitors of neutron stars in close binaries loose their 
hydrogen envelope as a result of mass transfer ( cf Section[2l) 
and th us end u p as typelb/c SNe. As show n bv lBrown et al.l 
(1 20011 ). see also lWellstein fc Langerl (|l999l ). stars which loose 
their envelope before or early during core helium burning (as 
in Case A/B RLO in a close binary) evolve as "naked" he- 
lium stars. There are four reasons why "naked" helium stars, 
at least at solar metallicity, develop small cores. Firstly, the 
lack of a hydrogen burning shell on top of the helium core 
prevents the core mass from growing during core helium 
burning. Secondly, the lack of this shell also prevents the 
convective part of the core from growing during core helium 
burning. In single stars, the growing convective core brings 
in fresh helium, which during late core helium burning is 
mostly used to convert carbon to oxygen. "Naked" helium 
stars therefore end core helium burning with a much higher 
carbon abundance (compared to single or wide orbit stars 
which evolve helium cores as "clothed"), and this will lead 
to lower mass iron cores as discussed above. Thirdly, helium 
stars emit strong winds as so called Wolf-Rayet stars, which 
reduces the final CO-core mass. And finally, not only are the 
core masses of "naked" helium stars reduced, compared to 
the situation of single or wide orbit stars, also the masses 
and the binding energies of the surrounding envelope are 
reduced, which leads to a corresponding reduction in the 
amount of material that falls back after the SN explosion 
IWooslev fc Weaver! [l99i l. All these effects lead to smaller 
neutron star masses in close binaries where the progenitor 
of the neutron star evolved through Case A/B RLO, com- 
pared to formation from a single star or wide binary (Case C 
RL O) evolution. 

iTimmes et al.l (| 19961 ) calculated remnant masses 
from explosio n s of the pre-SN helium star models of 
IWooslev et al.l (|l995l ). They estimate these models corre- 
spond to ZAMS masses up to 35 Mq for stars in close bi- 
naries which loose their hydr ogen envelope before or early 
during core helium burning. ITimmes et al. I l|l996t ) predict 
an upper limit for the gravitational mass of these neutron 
stars to be about 1.4 Mq (see their Fig. 6c for type lb SN). 
Depending on the details of the fall back of material these 
val ues could be as hi gh as ^ 1.6 M©. In fact, according 
to iBrown et al.l (1200 ll ) even stars initially as massive as 
40 — 60 Mq may end their lives as a neutron star if they loose 
their envelope early in their evolution (Case A or Case B 
RLO). 

At th e low end of progenitor masses, iPodsiadlowski et al.l 
l|2004h showed that electron capture supernovae from stars 
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in close binary systems may originate from a much broader 
mass range than in the single star (or wide binary) case due 
to the avoidance of the second dredee-up phase. 

In a recent paper lYoon et alj (|2010l ) investigated the 
evolution of stars from binaries leading to typelb/c SNe 
using the most recent theoretical models of (reduced) Wolf- 
Rayet winds. The result is larger final helium core masses 
which may lead to somewhat larger iron core masses and 
thus larger post-SN remnant masses. An investigation of this 
question is currently in progress. 

All in all, we see that the expected neutron star mass 
distribution from close binary stars is very different from 
that expected from single stars. In general, close binaries 
are thought to produce lower mass neutron stars. In this 
context, our conclusion that that the first mass transfer in 
the progenitor evolution of PSR J1614— 2230 was of Case C 
(a wide binary) is fully consistent with the high initial neu- 
tron star mass of 1.7 Mq we derived for this system above. 

4.3 Ramifications from PSR J1614 2230, 
Vela X-1 and the Blaclc- Widow pulsar 

Mass determinat ions of Vela X-1 l|Barziv et al.l l200ll : 
iRawls et all l201ll ) suggest that this neutron star has an 
observed mass of 1.77 ± 0.08 M©. The companion star to 
Vela X-1 is a BO. 5 lb supergiant (HD 77581) with a mass of 
about 23 Mq which implies that the present mass of the neu- 
tron star is very close to its birth mass. (Even a hypothetical 
strong wind accretion at the Eddington limit would not have 
resulted in accretion of more than about 10~^ Mq given the 
short lifetime of its massive companion). We therefore con- 
clude that not only was the neutron star in PSR J1614— 2230 
born massive (1.7±0.15 Mq) also the neutron star Vela X-1 
was born with a mass of at least 1.7 Mq. Furthermore, both 
of these neutron stars were produced from progenitors with 
a ZAMS mass above 20 Mq and these progenitors did not 
loose their envelope before core helium burning exhaustion 
(i.e. they evolved via Case C RLO prior to the SN - not 
to be confused with the RLO during the later X-ray phase 
discussed in Section[2]). Unlike PSR J1614-2230, which had 
a relatively low- or intermediate-mass companion. Vela X-1 
has a massive companion and this system is therefore not 
expected to have evolved through a CE prior t o the SN. 

A recent analysis bv Ivan Kerkwiik etlH (|201ll ) of the 
so-called Black- Widow pulsar yielded a mass of 2.4 Mq. Al- 
though the uncertainties of this result are rather large, such 
a mass would be difficult to explain if the neutron star was 
born with the canonical mass of about 1.4 Mq. 

4.4 Future observational constraints 

Future observations could push the empirical upper limit of 
the possible neutron star birth mass to even higher values. 
This could, for example, be achieved by measurements of a 
binary radio pulsar which reveal a massive neutron star in a 
double neutron star binary or, perhaps more likely, in a very 
tight binary with an 0-Ne-Mg WD. In both cases the evolu- 
tionary timescale of the progenitor of the last formed com- 
pact object would be so short that no substantial amount 
of mass could be accreted by the mildly recycled pulsar. 
Hence, in these cases the observed mass of the pulsar would 



be almost equal to its birth mass. Also detection of an even 
more massive neutron star than Vela X-1 in a HMXB system 
would push the limit upwards. 

We note that it is much more difficult and uncertain to 
estimate the birth mass of a pulsar in a close binary with 
a low-mass helium WD companion. Even in binaries with 
Porb < 1 day where the mass transfer rate is expected to 
have been sub-Eddington (e.g. Tauri s fc Savonii 3 1 19991 ) it 
is difficult posterior to estimate the effects on the accretion 
rate, and thus infer the birth mass of the neutron star, due to 
the poorly known occurrence of accretion disk instabilities 
and strength of the propeller effects - cf. further discussion 
in Paper II. 



5 CONCLUSIONS 

We have investigated the formation of PSR J1614— 2230 
by detailed modelling of the mass exchanging X-ray phase 
of the progenitor system. We have introduced a new ana- 
lytic parameterization for calculating the outcome of either 
a CE evolution or the highly super-Eddington isotropic re- 
emission mode, which depends only on the present observ- 
able mass ratio, q and the ratio between the initial and fi- 
nal donor star mass, k. Using a detailed stellar evolution 
code we calculated the outcome of a number of IMXBs 
undergoing Case A RLO. Based on the orbital dynamics 
and observational constraints on the stellar masses we find 
that PSR J1614— 2230 could have evolved from a neutron 
star with either a 2.2 — 2.6 Mq asymptotic giant donor star 
through a CE evolution (initiated by Case C RLO), or from a 
system with a 4.0 — 5.0 Mq donor star via Case A RLO. Sim- 
ple qualitative arguments on the evolution from the ZAMS 
to the X-ray phase suggest that Case A is the most likely of 
the two scenarios. The methods used in this paper, for RLO 
Cases A, B and C, could serve as a recipe for investigations 
of the progenitor system of other massive binary millisecond 
pulsars with heavy white dwarf companions to be discovered 
in the future. 

We conclude that the neutron star in PSR J1614— 2230 was 
born significantly more massive (1.7 ± 0.15 M©) than neu- 
tron stars found in previously known radio pulsar binaries 
- a fact which is important for understanding stellar evolu- 
tion of massive stars in binaries and the explosion physics of 
core collapse SNe. Finally, based on this high value for the 
neutron star birth mass we argue that the progenitor star of 
PSR J1614-2230 had a ZAMS mass of 20 - 25 Mq and did 
not loose its envelope before core helium exhaustion. 
In Paper II we continue the discussion of the formation of 
PSR J1614— 2230 in view of the spin-up process and include 
general aspects of accretion onto a neutron star during the 
recycling process and apply our results to other observed 
millisecond pulsars. 
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